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Abstract 

We propose a new non-perturbative method to compute derivatives of gauge 



> 

O ' coupling constants with respect to anisotropic lattice spacings (anisotropy 

■ coefficients), which are required in an evaluation of thermodynamic quantities 

from numerical simulations on the lattice. Our method is based on a precise 
measurement of the finite temperature deconfining transition curve in the 

Oh: 

lattice coupling parameter space extended to anisotropic lattices by applying 



the spectral density method. We test the method for the cases of SU (2) and 
SU(3) gauge theories at the deconfining transition point on lattices with the 
lattice size in the time direction Nt = 4 - 6. In both cases, there is a clear 
discrepancy between our results and perturbative values. A longstanding 
problem, when one uses the perturbative anisotropy coefficients, is a non- 
vanishing pressure gap at the deconfining transition point in the SU{3) gauge 
theory. Using our non-perturbative anisotropy coefficients, we find that this 
problem is completely resolved: we obtain Ap/T 4 = 0.001(15) and —0.003(17) 
on Nt = 4 and 6 lattices, respectively. 
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I. INTRODUCTION 



In order to study the nature of the quark-gluon plasma in heavy ion collisions and in the 
early Universe, it is important to evaluate the energy density e and the pressure p near the 
transition temperature of the deconfining phase transition. These quantities are defined by 
derivatives of the partition function in terms of the temperature T and the physical volume 
V of the system 

ldlnZ d\nZ 

V = T ^TT- (!) 



V dT- 1 ' * dV 

The lattice formulation of QCD provides us with a non-perturbative way to compute these 
quantities by numerical simulations. On a lattice with a size N% x N t , V and T are given by 
V = (N s a s ) 3 and T = l/(N t at), with a s and at the lattice spacings in spatial and temporal 
directions. Because N s and N t are discrete parameters, the partial differentiations in ([!]) are 
performed by varying a s and at independently on anisotropic lattices. 

The anisotropy on a lattice is realized by introducing different coupling parameters in 
temporal and spatial directions. For an SU (N c ) gauge theory, the standard plaquette action 
on an anisotropic lattice is given by 

s = -Ps E PiAx)-Pt Y, p d x )> ( 2 ) 

where P^ u (x) = ^-Re Tr{U IJ ,(x)U 1 ,(x + jl)!!^ (x+v)Ul(x)} is the plaquette in the (/i, v) plane. 
With this action, the energy density and pressure are given by fflJ 



e = —gr- I [a t — - 1 ((P s ) - (P> ) + I - £— ) ((P t ) - (P) ) \ , (3) 



P = ^ « P °) ~ ( P )o) + ^ ((Pt) ~ (P>o)} , (4) 

where (P s (t)) is the space(time)-like plaquette expectation value and (P) the plaquette 
expectation value on a zero-temperature lattice. Here, for later convenience, we have chosen 
a t and £ = a s /a t as independent variables to vary the lattice spacings, instead of a s and £ 
adopted in @. 



In order to compute e and p from eqs. (|3|) and (|4]) using numerical results from simulations, 
the values for the derivatives of gauge coupling constants with respect to the anisotropic 
lattice spacings 

d(3 a dfc d(3 s d(3 t 
at da~ t > at da~ t > flf' W {) 

which we call the anisotropy coefficients, are required. They can be computed from a 
requirement that, in the scaling region, the effects of anisotropy in the physical observables 
can be absorbed by a renormalization of the coupling parameters. Similar to the case of 
the renormalization group beta-function, the anisotropy coefficients do not depend on the 
temperature, because the renormalization is independent of the temperature. 

The calculation of these anisotropy coefficients in the lowest order perturbation theory is 
done by Karsch . However, the perturbative coefficients are known to lead to pathological 
results such as a negative pressure and a non-vanishing pressure gap at the deconfining 
transition in SU(3) gauge theory. Therefore, non-perturbative values of the anisotropy 
coefficients are required in order to study the thermodynamic quantities near the phase 
transition when N t is not sufficiently large. 

We are interested in the values of the anisotropy coefficients for isotropic lattices 
(/3 S = fit = (3, i.e. £ = 1) where most simulations are performed. In this case, we have 

( a *fff)£ =1 = ( a *lfr^ =1 = = ^ ca ^cnT' wnere a^^- is the beta-function at £ = 1, whose 
non-perturbative values are well studied both in SU{2) and SU(3) gauge theories Fur- 
thermore, a combination of the remaining two anisotropy coefficients is known to be related 
to the beta-function |2j by[] 

Therefore, only one additional input is required to determine the anisotropy coefficients for 
isotropic lattices. 



1 In H], a corresponding equation is given for {dj3 s n\ /0£)a s :fixed- 
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A non-perturbative determination of the anisotropy coefficients was attempted in 
Refs. [fH-|Tt| using a method that we call "the matching method" in the following. One 
first determines £ as a function of (3 S and fit by matching space-like and time-like Wil- 
son loops on anisotropic lattices, and then numerically determines dj/d£ at £ = 1, where 



7 = J/3 t /(3 s . Interpolation of the Wilson loop data at different sizes or interpolation of £ at 
different 7's using an Ansatz is required to evaluate dj/d^ at £ = 1. 

Alternatively, we can evaluate a non-perturbative value of pressure directly from the 
Monte Carlo data by "the integral method" [|ll|] : Assuming homogeneity expected when the 
spatial lattice size is sufficiently large, we obtain the relation p = —f, where / = —yhiZ 
is the free energy density, which can be evaluated by numerically integrating the plaquette 
difference (P s ) + (P t ) — 2(P)q in terms of (3 on isotropic lattices. The resulting value of the 
pressure, in turn, provides us with a non-perturbative estimate of an anisotropy coefficient 
@||. In actual numerical simulations, as the value of p in the confining phase and near 
the deconfining transition point is quite small compared with the magnitude of errors, it is 
difficult to determine the anisotropy coefficients near the transition point ||. 

In this paper, we propose a new method to directly compute the anisotropy coefficients 
at the deconfining transition point. Our method is described in Sec. |TJ. We test the method 
in the cases of SU (2) gauge theory in Sec. |T|. The more realistic case of SU(3) gauge theory 
is studied in Sec. |TV|. As an application of our non-perturbative anisotropy coefficients, we 
study the gaps for e and p at the SU(3) deconfining transition for N t = 4 and 6. A summary 
is given in Sec. [V|. 



II. METHOD 

Our method is based on an observation that, in the scaling region, the transition temper- 
ature T c = 1/ {N t a t (P s , Pt)} must be independent of the anisotropy of the lattice. Therefore, 
when we change the coupling constants along the transition curve in the (j3 s ,j3t) plane like 
(P s ,Pt) (Ps + d-Ps,Pt + d/? 4 ) on a lattice with fixed N t , the lattice spacing in the time 



4 



direction a t does not change: 

da t = ^d(3 s + ^d(3 t = 0. 
of3 s dp t 

We denote the slope of the transition curve at £ = 1 by r t \ 

(d(3 s \ 



(7) 



dp s _/fW\ / fdaA 
d(3 t = \d(5 t j i= J W s j i=l 



where we used an identity 



/ 



dPs dp t \ 

da t da t 

80s dpt 

V 9? at J 



<3g da t <9g dat 

d/3t df3 s dp s dfa 



on 



dat dat 

\ apt ap s J 



(8) 



dPs 



(9) 



Hence, the derivatives of (3 S and p t in terms of £ are expressed as 



3r t dP 



2(1 + r t ) da 
3 dp 



'OPA 3_ 

v 0£^ =1 2(1 + r t ) da 



(10) 



Introducing the conventional notation 7 = J p t / P s and /3 = \fP7Pt-, we obtain 



'dj s 



3 l-r t d/3 
a- 



' a, :fixed, e =l V ^ / a s :fixed, £=1 4/3 1 + T t da 

Finally, the customarily used forms for the anisotropy coefficients (Karsch coefficients 
are given by 



(11) 



'dg 



-2' 



a s :fixed, £=1 



H 2(1 + r t ) da 



2N, 



dg 



t 

w 



a s :fixed, £=1 
-2c — 1 „„A Q OAf „— 2 



2N r 



-P 



1 - 2r t dp 
a- 



(12) 



2(1 + r t ) da} 1 

where P s = 2N c g~' 2 ^~ 1 and Pt = 2N c gj~ z £ i . Therefore, when the value for the beta-function 
is available, we can determine these anisotropy coefficients by measuring r t from the finite 
temperature transition curve in the (P s ,p t ) plane.0 



2 A similar approach was proposed in [12]. 
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In order to determine the transition curve in the coupling parameter space, we compute 
the rotated Polyakov loop 

1 1 Nt 

L = z JpE W ^UU,(x,t) (13) 

S X c t = l 

as a function of (/3 s ,/3 t ), where z is a Z(N C ) phase factor (z Nc = 1) such that arg(L) G 
(— 7r/N c , 7r/N c ]. We define the transition point as the peak position of the susceptibility 
X = N*({L 2 } - (L) 2 ) in /3 for each fixed 7. 

We compute the coupling parameter dependence of \ m the (/3 s ,/3 t ) plane by apply- 



ing the spectral density method |T3[ extended to anisotropic lattices. This enables us to 
compute the anisotropy coefficients directly from simulations at £ ~ 1 without introducing 
an interpolation Ansatz. Another good feature of the spectral density method is that the 
method works well even with data obtained only on isotropic lattices. Therefore, we can use 
data from previous high statistic simulations performed on isotropic lattices, when the time 
histories of the Polyakov loop and space-like and time-like plaquettes are available near the 
transition point. 

Fitting the transition curve with a polynomial 

^max 

&(7) = E fn (7 - 1)" , (14) 
n=0 

with f n the fitting parameters, the slope r t is given by 

r _ ( mii) \ (d&/d 7 )e = i - Pc ( , 

1 U(&7)J e =i (d&/d7)*=i+&' 

where {df3 c /dq)^ = i = f\. The range of (3 and 7 in which the spectral density method is reliable 

is estimated by the condition that the statistical error for the reweighting factor (which is 

(e~ AS ) when the number of simulation points is one) is less than 0.5%. We confirm that 

the results are completely stable under a variation of n max when we restrict ourselves to the 

range discussed above. Choosing a range of 7 around 1 in such a way that the transition 

curve is almost straight, we use n max = 3 for the final results. 
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III. RESULTS FOR SU(2) 



We first test the method for the case of SU (2) gauge theory at the transition point /3 C for 
N t = 4 and 5. Although the method should work well with data only from isotropic lattices, 
in order to confirm it, we perform Monte Carlo simulations also on several anisotropic 
lattices for SU{2). On a 16 3 x 4 lattice, we perform simulations at (/3 S , (3 t ) = (2.300, 2.300), 
(2.302, 2.302), (2.296, 2.306), and (2.307, 2.298). On a 20 3 x 5 lattice, we simulate at 
(P 8 ,Pt) = (2.373,2.373), (2.375,2.375), (2.380,2.370) and (2.368,2.378). At each on 
the N t = 4 (5) lattice, we accumulate 500,000 (1,250,000) configurations, each separated 
by 10 heat-bath sweeps, after thermalization. The statistical errors are estimated using the 
jackknife method with the bin size of 1000 configurations. We confirm that the errors are 
stable under a wide variation of the bin size around this value. 

Computing the susceptibility in the (/3 s ,/3 t ) plane using data at each simulation point, 
we check that the results agree well with each other, i.e. the results for the susceptibility 
from isotropic lattices coincide with the results from anisotropic lattices. For the rest of this 
section, we combine the results for all four (/3 S , (3 t ) combinations to compute the susceptibility 
with the spectral density method. In Fig.0, we plot the susceptibility for N t = 4 at 7 = 0.995, 
1.000, and 1.005. The results for the peak position (3 C of the susceptibility computed at 
various values of 7 are summarized in Fig. |2] for N t = 4 and 5. 

Fitting the results for the transition curve, we obtain the values for (5 C and rt at £ = 1, 
as summarized in Table [j| Combining the values of r t with a result of the SU(2) beta- 
function H at /3 C (£ = 1), we obtain the anisotropy coefficients (|i~l~D and (|12|) . The results 
are summarized in Table 0. Because no errors for the beta-function are given in ||, we 
disregard their contribution to the errors of the anisotropy coefficients. 

In Fig. |3|, we compare our results for the Karsch coefficients with the results of the 
perturbation theory (dot-dashed curves) and the integral method (dotted curves) [§]]. We 
find significant discrepancies between our results and the results of the perturbation theory. 
On the other hand, our results are consistent with the results from the integral method. 
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IV. RESULTS FOR SU(3) 



Let us now study the more realistic case of the SU(3) gauge theory. We analyze the 
high statistic data for the SU(3) gauge theory obtained by the QCDPAX Collaboration 
[[14]]. Simulations were performed at the deconfming transition point for N t = 4 and 6. For 
N t = 4, the lattice sizes are 24 2 x 36 x 4 and 12 3 x 24 x 4, with 712 000 and 910 000 pseudo 
heat-bath iterations, respectively. For Nt = 6, data on 36 2 x 48 x 6, 24 3 x 6, and 20 3 x 6 
lattices with 1 112 000, 480 000, and 376 000 iterations are available. The Polyakov loop and 
the plaquettes are measured every iteration. Details of the simulation parameters are given 



in [|14l . For the bin size in the jack-knife analysis, we adopt the same values as in |14 



A. Anisotropy coefficients 

The results for the susceptibility on the largest spatial lattices are given in Figs. [| and ||. 
Because the transition is of first order for SU(3), the peak of the susceptibility is quite clear 
when the spatial lattice size is large enough, as shown in Figs. ^ and |5|. (Note the difference 
in the vertical scales between Figs. [1] and 

Our results for the slope r t are summarized in Table III]. Except for the case of the 24 3 x 6 
lattice where the simulation point is slightly off the transition point, the errors become larger 
with decreasing spatial volume, because the peak of the susceptibility becomes less clear on 
small lattices. From Table p| , we find that the slopes at Nt = 4 with different spatial 
lattice volumes completely agree with each other. As shown in Figs. |] and [|, the peak of 
the susceptibility for Nt = 6 is less sharp compared with that for N t = 4 with the same 



relative spatial volume (N s /N t ) due to the fact that the transition is weaker for A f = 6 []14 
Therefore, with comparable statistics, r t has a larger statistical error for Nt = 6. Unlike in 
the case of N t = 4, the central values for the slope for N t = 6 given in Table |TTT| vary with the 
spatial volume by about one standard deviation. However, because the volume dependence 
is not uniform, we consider that it is caused by statistical fluctuations. We use the values 



S 



obtained on the largest spatial lattices for our final results. 

Our results for the anisotropy coefficients are summarized in Table [IV]. For our final 
results, we adopt the beta-function computed from a recent string tension data by the SCRI 
group 0. See a subsection below for a discussion about the influence on the results from 
the choice of the beta-function. 



B. Pressure gap and latent heat 

As an application of our non-perturbative anisotropy coefficients, we reanalyze the ther- 
modynamic quantities e and p at the deconfining transition point using the plaquette data 



by the QCDPAX Collaboration [[L4j]. In terms of the slope r t and the beta-function, the 



conventional combinations e — 3p and e + p are given by 



(e - 3p)/T 4 = -3A?a^ {(P s ) + (P t ) - 2(P) }, (16) 



(e + p)/T 4 = 3A? ^ {(P.) - (P t )}. (17) 

At a first order transition point, we have a finite gap for energy density, the latent heat, but 
expect no gap for pressure. It is known that the perturbative anisotropy coefficients have 
a difficulty which leads to a non- vanishing pressure gap at the deconfining transition point: 
Ap/T 4 = -0.32(3) and -0.14(2) at N t = 4 and 6 @ . 



New values for the gaps in e and p using our non-perturbative anisotropy coefficients are 
summarized in Table |VJ. For the pressure gap, we obtain 

0.001(15) foriV t = 4, 
Ap/T 4 = (18) 

I -0.003(17) for N t = 6. 

We find that the problem of non-zero pressure gap is completely resolved with our non- 
perturbative anisotropy coefficients. 
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C. Choice of the beta-function 



In Table we study the influence of the choice of the beta-function on the anisotropy 
coefficients. We compare (i) the beta-function computed from a recent string tension data 
by the SCRI group §, (ii) that from a MCRG study by the QCDTARO Collaboration §, 
and (iii) that from a study of (3 c {N t ) by the Bielefeld group ||. The SCRI beta-function is 
computed using a fit of the string tension for 5.6 < (3 < 6.5. We note that the QCDTARO 
beta-function is based on a fit of mean-field improved gauge coupling constant using the 
results of plaquette at (3 > 5.8; i.e. f3 c (N t = 4) « 5.69 is slightly off the range of validity P,|T5|. 
Also the beta- function by the Bielefeld group seems to be problematic around /3 c (N t = 4), 
because it is largely affected by the data of (3 c (N t = 3) where we cannot expect universal 
scaling. Accordingly, the beta-function of the Bielefeld group shows a systematic deviation 
from the data of a MCRG study at (3 < 6 §. 

These beta-functions are plotted in Fig. ||. At fl c (N t = 6), different beta-functions coincide 
with each other within 5%, while, at (3 c (N t = 4), they vary by about 20%. Because only 
the SCRI beta-function is reliable at P c (N t = A) as discussed in the previous paragraph, we 
adopt the SCRI beta-function for our final results. 

In order to compare the anisotropy coefficients from different references, however, it is 
important to check the effect of the beta-function on the results. From Table [TV], we see 
that the results for the anisotropy coefficients using different beta-functions agree well with 
each other at N t = 6. At N t = 4, however, the anisotropy coefficients depend very much 
on the choice of the beta-function. Accordingly, we find that the results for the latent heat 
are consistent with each other at N t = 6: Ae/T 4 = 1.569(40), 1.539(39), and 1.515(38) 
with SCRI, QCDTARO, and Bielefeld beta-functions, respectively. At Nt = 4, we find a 
sizable dependence on the choice of the beta-function: Ae/T 4 = 2.074(34), 1.877(30), and 
2.265(37) using SCRI, QCDTARO, and Bielefeld beta-functions. For the pressure gap, on 
the other hand, because the beta-function appears only as a common overall factor in ( |I~6"|) 
and (|T7|), the conclusion that Ap vanishes with our anisotropy coefficients does not depend 
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on the choice of the beta-function. 



D. Comparison with other methods 

In Fig. 0, we summarize our results for the Karsch coefficients together with previous 
values; the perturbative results 0, results from the integral method ||, and those from 
the matching of Wilson loops on anisotropic lattices PJlQl. No errors are published for the 
results from the integral method. We find that all non-perturbative methods give values 
which deviate from the results in the perturbation theory. 

Comparing the results from different non-perturbative methods, we find that, although 
the deviations from the perturbation theory are roughly consistent with each other, the cen- 
tral values are different by more than three standard deviations, when we take the published 
errors. 

We think that one origin of the variation among different methods at (3 c (N t = 4) is the 
beta-function. Note that the results from Refs. || (matching method) and |J (integral 
method) are computed using the beta-function of the Bielefeld group, while our results 
and the results from Ref. (matching method) are using the SCRI beta-function. From 



Table [TV], we note that, if we adopt the beta-function of the Bielefeld group, our results are 
consistent with those of Ref. || at (3 c (N t = A). 

At (3 c (N t = 6), on the other hand, the difference in the results is not due to the beta- 
function, because the systematic error due to the choice of the beta-function is small as 
discussed in the previous subsection. In order to see this, we study d'j/d!;, which can be 
computed without using the beta-function in the matching method. The values of d'y/dl; 
obtained in Ref. [|1(| are reported to be consistent with those from the integral method 0, 



but are different to another result from the matching method ||. Performing a quadratic 
interpolation in (3, we find dj/d£ ~ 0.64(1) (lj, 0.66(2) §, and 0.74(2) § at (3 c (N t = 6). 



Our result 0.707(10) given in Table [IV] is around the center of these values. A careful study 
of systematic errors in each method is required to understand the variation between different 
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methods. 



V. CONCLUSIONS 

We have computed the anisotropy coefficients for the 577(2) and SU(3) gauge theories 
by measuring the transition curve of the deconfining transition in the (f3 s ,f3t) plane. One of 
the essential ingredients of our approach is the application of the spectral density method, 
that enables us to determine the anisotropy coefficients directly from simulations at £ ~ 1. 
We note that the spectral density method is useful to avoid interpolation Ansatze also in 
the matching method. 

Our non-perturbative results for the anisotropy coefficients are summarized in Tables [II 
and |V[ Our results shown in Fig. |7| suggest that the Karsch coefficients converge to the 
perturbative values slightly faster than that suggested by the central values from Refs. || 
and [[KJ. Applying the results for SU(3), we reanalyzed the thermodynamic quantities at 
the deconfining transition point on N t = 4 and 6 lattices. We obtain vanishing pressure gaps 
with our non-perturbative anisotropy coefficients, thereby solving a longstanding problem 
of non-zero pressure gap with the perturbative coefficients. 

We are grateful to O. Miyamura, A. Nakamura and H. Matsufuru for useful discussions 
and sending us the data for the QCDTARO beta-function. We also thank A. Ukawa, T. 
Yoshie, Y. Aoki, T. Kaneko, R. Burkhalter and H.P. Shanahan for helpful suggestions and 
comments. This work is in part supported by the Grants-in-Aid of Ministry of Education, 
Science and Culture (Nos. 08NP0101 and 09304029). SE is supported by the Japan Society 
for the Promotion of Science. 
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TABLES 

TABLE I. Results for (5 C and the slopes at £ = 1 in the SU(2) gauge theory. The column 
"7-range" is for the range of 7 used in the fit for the slope d(5 c /d , y. 



lattice 


(3c 


7-range 


df3 c /d-f 


n 


16 3 x 4 
20 3 x 5 


2.30177(9) 
2.37430(8) 


0.995 - 1.005 
0.995 - 1.005 


-0.370(12) 
-0.312(15) 


-1.384(14) 
-1.303(17) 


TABLE II. SU (2) anisotropy coefficients at £ = 1 using 
by the Bielefeld group Q. 


the beta-function adg~ 


~ 2 /da obtained 


lattice 


d^/di 


c s 




adg~ 2 /da 


16 3 x 4 
20 3 x 5 


0.683(21) 
0.725(35) 


0.203(12) 
0.182(21) 


-0.161(12) 
-0.144(21) 


-0.08439 
-0.07544 


TABLE III. The 


same as Table [ 


| for SU (3) using the data 


by the QCDPAX Collaboration [14|. 


lattice 


0c 


7-range 


dpc/dj 


n 


24 2 x 36 x 4 
12 2 x 24 x 4 


5.69245(23) 
5.69149(42) 


0.9975 - 1.0025 
0.995 - 1.005 


-0.5193(23) 
-0.5183(52) 


-1.2008(10) 
-1.2004(22) 


36 2 x 48 x 6 
24 3 x 6 
20 3 x 6 


5.89379(34) 
5.89292(87) 
5.8924(14) 


0.999 - 1.001 
0.999 - 1.001 
0.9975 - 1.0025 


-0.5844(83) 
-0.542(33) 
-0.622(34) 


-1.2201(35) 
-1.202(14) 
-1.236(14) 
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TABLE IV. SU(3) anisotropy coefficients at £ = 1, using the values for the beta-function 
adg~ 2 /da by the SCRI group §, the QCDTARO Collaboration g, and the Bielefeld group §. 
For our final results, we take the values obtained on the largest spatial lattices using the SCRI 
beta-function. Because the errors for the beta-function are not given in the papers, we disregard 
their contribution to the errors of the anisotropy coefficients in this table. See text for details. 
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U.OoZZ v ZO J 


— U.O400 v ZO ) 


— U.U 1 lUo 


ovarii 




U.00 i v ZO J 


u.^ooy v zo j 


— U.^UO / v ZO J 
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U.O / Zo v OU J 
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fl 9Q1 
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— u.o^o^oy j 


— u.u / uy i 






U.00 i O^Ou J 


U.^oOU v OO J 


— u.^uoy x oo j 


— U.UU^ilO 


V / V 1 / 1 . \ I V V / 




U.O 1 OO^Do J 


U.OZoo v 04 ) 


— u.zyuu v 04 ) 


— U.U 1 / 01 


Dieieieia 


36 2 x 48 x 6 


0.7068(100) 


0.3109(98) 


-0.2650(98) 


-0.09179 


SCRI 




0.6936(98) 


0.3235(96) 


-0.2784(96) 


-0.09008 


QCDTARO 




0.6826(96) 


0.3340(95) 


-0.2897(95) 


-0.08864 


Bielefeld 


24 3 x 6 


0.762(47) 


0.257(46) 


-0.211(46) 


-0.09172 


SCRI 




0.747(46) 


0.271(45) 


-0.226(45) 


-0.08999 


QCDTARO 




0.736(45) 


0.282(45) 


-0.237(45) 


-0.08857 


Bielefeld 


20 3 x 6 


0.663(36) 


0.354(35) 


-0.308(35) 


-0.09167 


SCRI 




0.651(35) 


0.366(35) 


-0.321(35) 


-0.08994 


QCDTARO 




0.640(35) 


0.375(34) 


-0.331(34) 


-0.08853 


Bielefeld 
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TABLE V. Gaps for thermodynamic quantities in the SU(3) gauge theory at the deconfining 
transition point using our non-perturbative anisotropy coefficients. Plaquette data are taken from 
Ref. [|l4|]. The low temperature hadronic phase (had) and the high temperature quark-gluon-plasma 
phase (QGP) are separated as described in ]i~4|| . We reanalyze (e — 3p)/T 4 also, using the SCRI 
beta-function. 



lattice 


24 2 x 36 x 4 


36 2 x 48 x 6 


P 


5.6925 


5.8936 


A(e + p)/T 4 


2.075(42) 


1.565(51) 


A(e- 3p)/T 4 


2.072(43) 


1.578(42) 


Ae/T 4 


2.074(34) 


1.569(40) 


Ap/T 4 


0.001(15) 


-0.003(17) 
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FIG. 1. Polyakov loop susceptibility in the SU(2) gauge theory on a 16 3 x 4 lattice at 7 = 0.995, 
1.0 and 1.005. Errors are estimated by a jackknife method. 
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(a) 




2.365 2.37 2.375 2.38 




FIG. 2. Polyakov loop susceptibility in the SU{2) gauge theory as a function of {(3 s ,(3t) 
obtained on (a) 16 3 x 4 and (b) 20 3 x 5 lattices. Simulation points are shown by filled circles. 
The bold lines represent the peak position of the susceptibility and the dashed lines their errors. 
The magnitude of the susceptibility is shown by tone for the range (a) 8.2 < x < 10-4 and (b) 
9.0 < x < 11-2, respectively, where different tone corresponds to a difference Ax = 0.2. 
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FIG. 3. Anisotropy coefficients c s and ct for the SU{2) gauge theory. Our non-perturbative 
results are given by filled circles. The dot-dashed curves are the results of the perturbation theory 
Qj. The dotted curves are the results from the integral method m. No errors are published for 
these curves. 
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(a) 





FIG. 4. Polyakov loop susceptibility in the SU(3) gauge theory obtained (a) on the 24 2 x 36 x 4 
lattice at 7 = 0.9975, 1.0 and 1.0025, and (b) on the 36 2 x 48 x 6 lattice at 7 = 0.999, 1.0 and 
1.001. 
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FIG. 5. The same as Fig. |2] for the SU(3) gauge theory on (a) 24 3 x 36 x 4 and (b) 36 2 x 48 x 6 
lattices. The range of x plotted and the width Ax for a tone are (a) 0.0 - 45.0, 5.0 and (b) 2.5 - 
22.5, 2.5, respectively. 
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FIG. 6. Non-perturbative beta-functions in the SU(3) gauge theory. 
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FIG. 7. Anisotropy coefficients in the SU{3) gauge theory. Our non-perturbative results are 
given by filled circles. The dot-dashed curves are the results of the perturbation theory The 
open squares are those from a matching of Wilson loops |J. Open triangles and thin lines are the 
results of a matching method combined with the SCRI beta-function [||. The dotted curves are 
the results from the integral method (BJ. No errors are published for the results from the integral 
method. 
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